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Mh' Abstract I provide a broad overview of the basic theoretical paradigms of black hole 

■^L . accretion flows. Models that make contact with observa tions continue to be mostl y 

based on the four decade old alpha stress prescription of IShakura &: SunvaevI (| 19731 1 . 

and I discuss the properties of both radiatively efficient and inefficient models, including 

their local properties, their expected stability to secular perturbations, and how they 

might be tied together in global flow geometries. The alpha stress is a prescription 

M— I . for turbulence, for which the only existing plausible candidate is that which develops 

from the magnetorotational instability (MRf). I therefore also review what is currently 

_7 , known about the local properties of such turbulence, and the physical issues that have 

f-«) ' been elucidated and that remain uncertain that are relevant for the various alpha-based 

I black hole accretion flow models. 
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1 Introduction 

>' 

0^ ' Accretion is the very process that allows black hole sources to emit electromagnetic 

[^^ , radiation and other forms of energy. Because black holes are so small in size compared 

to the spatial scale of their sources of fueling, and because centrifugal forces on matter 
of given angular momentum increase more rapidly (oc R~ ) than gravity (oc _R~ ) as 
nJ , one moves inward in radius 7?, accretion is generally believed to be a process involving 

>-^ ' rotationally supported flows. Matter in such a flow must lose angular momentum in 

order to move inward and release gravitational binding energy. It is the nature of the 
angular momentum loss mechanism, and the process whereby gravitational binding 
energy is converted into observable forms of energy, that are the two central questions of 
black hole accretion theory. At least three mechanisms have been proposed for angular 
momentum extraction: 
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(1) External stresses associated with larg e scale magnetic fields in a magnetohydrody 



namic (MHD) outflow. This mechanism i Blandford &: Pavng|l982l l may be relevant in 



low luminosity sources where accretion power may be largely converted into mechanical 
power in outflows. It may also be relevant in resolving the f ueling and self-gr avity prob- 



lems in the outer accretion flows in active galactic nuclei ( GoodmanI 120031 ) . Whether 



and how large scale magnetic fields can be created remains an open question, however. 
(2) Magnetorotational (MRI) turbulence. Such turbulence is generic for plasmas that 



are sufficiently electrically conducting and not too strongly m agnetized i Balbus fc HawlevI 
ll99l':'Ha wlev fc Balbuslll99ll :lBalbus fc Hawlevlll992l .ll998Y Because turbulence is in- 
herently dissipative, this process is almost certainly relevant for sources whose power 
output is dominated by thermal radiative emission mechanisms. 

(3) Nonaxisymmetric waves and shocks. Nonaxisymmetric (e.g. spiral) waves can trans- 
port angular momentum outward through the flow. Such waves can also transport en- 
ergy away from the region where gravitational binding energy is released, depositing 
it elsewhere. Waves are almost certainly relevant in disks around supermassive black 
hole binaries, and also in the outer, self-gravitating parts of disks in active galactic 
nuclei. They probably also play a role in the outer parts of black hole X-ray binary 
disks due to tidal excitation by the companion star. Nonaxisymmetric shocks can also 
play an important role in the inner regions o f accretion flows whos e angular momenta 
are misaligned with the black hole spin axis ( Fragile fc BlaedbOOa ). 



Among these options, only the second - MRI turb ulence - is a mechanism tha t might be 
describable by th e classical alpha prescriptio n of IShakura fc SunvaevI (119731 ). at least 



in some aspects ( Balbus fc Papaloizoulll999l ). The angular momentum transporting 



stress wji^ in the turbulence is given by local space and time averages of correlated 
fluctuations in radial (R) and azimuthal (0) fluctuations of velocity v (the Reynolds 
stress) and magnetic fleld B (the Maxwell stress), 

t^B4 = ( P'"BSy4, j^ ) , (1) 

where p here is the mass density and Sv^ is the local deviation of the azimuthal veloc- 
ity component from the mean background shear flow. Maxwell stresses are generally 
larger in magnitude than the Reynolds stresses by factors of at least several. I say that 
the total stress might be describable by the classical alpha prescription because these 
stresses appear to be mostly local in the sense that simulations show that radial corre- 
lations in stress drop rapidly on scales larger than the local disk scale height. However, 
as I discuss in section 3.1 below, there remain correlations on larger radial scales. 

This article provides a broad overview of alpha-based models of black hole accretion 
flows, focusing on structure, dynamics, and thermodynamics. These models continue to 
dominate theoretical efforts to explain observations, but a slow revolution is occurring 
as simulations of MRI turbulence, both local and global, continue to become more 
powerful and to incorporate more and more of the relevant physics. This article will 
also discuss what has been learned recently from local, shearing box simulations of MRI 
turbulence as this pertains directly to some of the alpha-based modeling. A review 
of global simulations can be found in Chapter 2.4. Spectral modeling of accretion 
flows is discussed in Chapter 2.3. I will also mainly focus on accretion rather than the 



formation of jets and outflows here, though jets and outflows are clearly important 
(both observationally and theoretically, in certain flow states). See Chapter 5.3 on jet 
launching mechanisms. 



2 Hydrodynamic disk models with the alpha prescription 

Decades of theory and models of black hole accreti on flows have criticall y relie d on the 
alpha prescription for a local stress introduced by IShakura fc SunvaevI ( 19731 1. There 



are numerous variants of this prescription which produce order unity changes in the 
definition of a, and one of the most common is 



w 



B.4> 



= aP, (2) 



where P is the thermal pressure. Most models have assumed that this is the total 
thermal pressure (gas plus radiation), but prescript ions in which the stress is taken to 
be proportional to just the gas pressure alone (e.g. ISakirnot o fc Co roniti 1981) or the 
geometric mean of the gas and total thermal pressures fe.g. lTaam fc Linlll984l ') have also 



been suggested. However, as illustrated in Figure[TJ recent radiation MHD simulations 
of MRI turbulence find that the stress scales best with total thermal pressure, at least 
on long time scales (lOhsuea et al.ll2009l : iHirose. Blaes fc Krolikll2009h . 

The alpha prescription ^ is usually used to solve for the radial structure of 
vertically-integrated geometrically thin or slim accretion disks, in which case it en- 
ters the equations through the vertically-integrated stress: 



Wr^ = / WR^dz ~ 2HaP, (3) 

<y — oo 

where P is now some vertically averaged thermal pressure, of order the midplane pres- 
sure, and H is the vertical half-thickness of the disk. This is consistent with MRI 
simulations, but the prescription is also occasionally used even more locally by assum- 
ing that the vertical profiles of stress and dissipation at a given radius are proportional 
to the local vertical profile of thermal pressure. As we discuss below in section 3.2, 
this is not consistent with vertically stratified simulations of MRI turbulence, which 
generally have vertical profiles of stress that are broader than the thermal pressure 
profile. Alpha defined locally would therefore increase rapidly outward from the disk 
midplane. 

2.1 Local thermal equilibria and secular instabilities 

Virtually all (non-simulation-based) models of black hole accretion flows are based on 
vertically integrated hydrodynamic equations. These models often neglect the possi- 
bility of significant losses of mass, angular momentum, and energy in outflows and 
jets, though some models do attempt to include them with various prescriptions, par- 
ticularly in advection-dominated flows which we will come to shortly. As discussed in 
Chapter 5.1, neglect of outflows is likely to be a bad approximation in some sources 
and accretion states. Nevertheless, if we adopt this assumption for simplicity, then for 
stationary flows, the conservation laws of mass, radial momentum, angular momentum, 
and internal energy can be written as 

M = 2ttREv, (4) 
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Fig. 1 Time-averaged values of the ratio of spatially averaged stress to various measures of 
spatially averaged thermal pressure, i.e. the lShakura &: SunvaevI ||1973| ) parameter alpha, as a 
function of the time-averaged ratio of spatially averaged radiation pressure to spatially averaged 
gas pressure, for a number of radiation MHD, vertically stratified, shearing box simulations of 
MRI turbulence. Black, green and blue points are the results for thermal pressures defined to 
be the total (radiation plus gas) pressure, the geometric mean of the total and gas pressures, 
and the gas pressure alone, respectively. Horizontal and vertical error bars on all points indicate 
one standard deviation in the respective time-averages. The horizontal black line is the average 
alpha value of the total pressure prescription (black) points, while the green and blue curves are 
what would result if the total pressure prescription were correct, but one nevertheless insisted 
on defining alpha in terms of the other thermal pressure definitions used in the green and blue 
points, respectively. The stress prescription that is most consistent with the simulation data 
is one in which the total thermal pressure is used, though it is perhaps noteworthy that the 
alpha values in the gas pressure dominated simulations are consistentl y higher t han the alpha 
values in the radiation pressure dominated simulations. (Updated from lHirose. Blacs fc KrolikI 
120091 .') 
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where we have neglected general relativity for the purposes of physical transparency. 
Here p ~ S/{2H) is a vertically-averaged density, S is the surface mass density, O is the 
fluid angular velocity which may differ from the test particle (Kepler) angular velocity 
i^K) ^ ~ ^R is the fluid specific angular momentum, v is the inward radial drift speed, 
[/ is a vertical average of the internal energy per unit mass, Q~ is the radiative cooling 



rate per unit surface area on each face of the disk, Q^ = —{l/2)Wji^RdQ/dR is half 
the turbulent dissipation rate per unit surface area, and Qadv is half the inward radial 
advection of heat per unit surface area. Vertical hydrostatic equilibrium implies that 
the vertical half-thickness of the disk is H ^ (P/p) /^K- 

Once one adopts the alpha prescription ((2]), together with an equation of state and 
opacities and/or optically thin cooling functions, it is possible to solve these equations 
with assumed boundary conditions to derive the radial profiles of vertically-averaged 
fluid variables in the flow. Such models generally invoke a regularity condition at an 
inner sonic point and/or a no-torque inner boundary condition at, for example, the 
innermost stable circular orbit, although magnetohydrodynamic stresses can be im- 
porta nt once one enters the plun ging region near th e bl ack hole (Ganmiie 1999: Kr ^oliM 
Il999l ). (See Chapter 2.4 and, e.g.. iPenna eTaLlboiol and lNoble. Krolik fc Hawlev II2OI0I 



for recent simulation work on this issue for geometrically thin disks.) Another approach 
is to consider radii much larger than the gravitational radius _Rg = GM/c and invoke 
self- similarity by assuming a constant ratio of advective cooling over turbulent dissi- 



pation Qadv/Q"*" (|Naravan fc Yilll994l '). 



For a fixed black hole mass, models that are stationary generally depend on a 
number of chosen parameters, the most important being the accretion rate M which 
is everywhere constant through the flow (remember, we are neglecting outflows here). 
A number of possible equilibria have been discovered in this way, and the primary 
method of choosing which ones are physically realizable in nature has been to check 
if they are stable to secular perturbations. The growth rates of such instabilities are 
related to one of two characteristic time-scales. The flrst is the thermal time, defined 
as the characteristic heating time 

Thermal instabilities, in which a local patch of the flow undergoes runaway heating or 
cooling, generally grow on this time scale. The second time scale is the inflow time, i.e. 
the time it would take for turbulent stresses to cause a fluid element to drift inward 
over a distance comparable to its current radius. 



R EQR^ n / R 



2Hwr^ afi^ 



(5)- 



^^ Viscous" or inflow instabilities tend to grow on this time scale, where I have enclosed 
the term "viscous" in quotation marks here (only) to emphasize that it is turbulent 
stresses, not microscopic viscosity, that play a role here. In geometrically thin accretion 
disks, where dynamical equilibrium on the inflow time requires the angular velocity and 
specific angular momentum of fiuid orbits to be a function of radius and not of time, 
one can write down a tim e-dependent diffusion equation for the evolution of the sur face 



one can write down a tim e-dependent dinusion equation tor the evolution 01 the sur t 
mass density in the flow (JLvnden-Befl fc Pringijl 19741 : iLightman fc Eardlevlfl973 '). 

(10) 
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where £' = {d/dR){R SI) is the radial derivative of specific angular momentum. In 
this case, instabilities happen if Wf{,p is inversely related to the surface mass density, 
as this equation then corresponds to a diffusion equation with a negative diffusion 
coefficient. Perturbations in surface mass density would then tend to grow, rather than 
be smoothed out, by this anii-diffusion. 



FoUowin K theoretical work on dwarf nova outbursts in cataclysmic variables (e.g. 
ISmakI 1 19841 ). it has proved convenient to depict local thermal equilibrium (Q = 
Q~ + Qadv) solutions at one p articular radius in a diagram of accretion rate vs. local 
surface mass density. Figure El ( Chen et al.lll995l ) depicts the topology of the space of 



such solutions at a radius of lOGM/c around a ten solar mass Schwarzschild black 
hole. Each curve corresponds to a different chosen value of a, as labelled. The locations 
of these curves in this graph can shift considerably, depending on the particular physics 
being included in the models, but the topological structure is robust. Naively, we expect 
equilibrium curves with negative slopes in this diagram to be viscously unstable, as the 
vertically integrated stress M^^^ is proportional to the dissipation rate per unit area Q 
which in turn is proportional to the accretion rate Al , because, after all, it is the inflow 
of matter which is the source of accretion power. The negative sloped curve portions in 
the right of this diagram near Eddington accretion rates are radi ation dominated, geo- 



metr ically thin disks, and are viscously u nstable by this criterion (iLightman fc Eardlev 
1974). They are also thermally unstable ( Shibazaki fc Hoshilll975l : IShakura fc Sunvaev 



19761 ) ■ as can be seen because Q exceeds (is less than) Q~ + Qadv above (below) 
these curve portions. Hence a perturbation upward (downward) from this curve will 
cause runaway heating (cooling), moving it away from the equilibrium curve. Simi- 
larly, the middle bold line in the lower portion of this diagram, which corresponds to 
a hot, optically thin flow where turbulent dissipation is balanced by radiative cool- 
ing is visco usly stable but thermally unstable. This solution was first discovered by 
Shapiro. Li ghtman & Eardlcy (1976). 

This leaves three regions of the diagram which appear to correspond to thermally 
and viscously stable solutions. The lower right set of curves correspond to the gas 
pressur e dominated regime of the o riginal geometrically thin, optically thick accretion 



disks of lShakura fc SunvaevI j 19731 ). At higher accretion rates radiation pressure starts 



to dominate these geometrically thin solutions and the curves bend over to the un- 
stable negative slopes. At still higher accretion rates, the inflow time becomes shorter 
than the cooling time, so that heat is advected inward. This advection is stabilizing 
l Abramowicz et al.lll988l ) and, because the scale height of the disk can become quite 



large, these flows have been dubbed "slim disks". (Note that for larger radii, the op- 
tically thick, geometrically thin disk curves would also bend toward unstable negative 
slopes and then back toward positive stable slopes as one lowers the accretion rate and 
passes through the transition where ionized hydrogen becomes neutral. This unstable 
branch is responsible for the transient outburst behavior observed in many black hole 
and neutron star X-ray binaries, as well as dwarf novae in accreting white dwarf sys- 
tems. We will have little to say about this instability here, as it generally occurs in the 
outer, less luminous portion of the disk in accreting black hole systems, but i t is crucial 
for explaining the phenomenology of black hole transients. See lLasotall200ll for a good 
review of the theory. Note, however, that these outbursts provide some of the only 
constraints on the levels of turbulent stress in accretion disks, with q ~ 0.1 — 0.3 in the 
outburst phase. This is significantly higher than that measured in local shearing box 
simulations with no net vertical magnetic flux, as illustrated in Figure [l] above. It may 
be that the character of the turbulence changes when one is so close to the regime of 
hydrogen ionization, or it may be that external magnetic fl ux is necessary to explain the 
observations. This is a significant unsolved problem, e.g. iKing. Pringle fc Livid 120071 : 
iKotkofc Lasotall2012l .) 

In addition to these two optically thick accretion disk solutions, a third set of stable 
solutions exists which is optically thin and which only exists at low accretion rates. 
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Fig. 2 Thermal equilibrium curves of particular accretion flow models around a ten solar 
mass black hole at a particular radius lOGM/c^ on the accretion rate (scaled wi th Mei the 
Eddington luminosity divided by c^) vs. surfac e mass density (17) plane, from [Chen et all 
Ill995l ) . Each curve is labelled by the value of the IShakura &: SunvaevI lll973f) stress parameter 
a chosen in the model. The locations of these curves in this plane can change considerably 
depending on the physics being incorporated and how it is treated in these models - see 
IChen et al.l lll995l ) for more details. The topology of the curves is, however, reasonably robust. 



provided a is not too high. Here the cooling time significantly exceeds the infall time, so 
that advection again balances turbulent dissipation, and this apparently also stabilizes 



the flow (llchimaruliyilNaravan fc Yill994l : lAbramowicz et al.ll 19951 : IChen et al.lll995l : 



iNaravan fc Yilll995alll995bl '). Such flows are known as Advection Dominated Accretion 
Flows (ADAF's, a term which can also be applied to the optically thick, radiation 
dominated slim disks) or radiatively inefficient accretion flows (RIAF's). 

While the above arguments suggest that these three sets of solutions are thermally 
and viscously stable, attempts to rigorously demonstrate this involve subtle issues, 
particularly in the two advection dominated solutions, and the situation is in fact 
not entirely clear. Advection dominated solutions are fairly geometrically thick {H ~ 
R), so the thermal and inflow time scales are comparable, and there is no longer 
a clean separation of thermal and viscous instabilities. Thermal pressure is also not 
dynamically negligible, so that thermal perturbations can alter the speciflc angular 



momentum distribution and the surface mass density in the disk, even without mass 
diffu sion due to turbulent stresses. A t tempts to tackle this prob l em have been made 
fe.g iKato. Abramowicz fc Chenllliii : IWull 19961 : iKato et aL|[l997l : IWull 19971 : lYamasakil 
119971 ). but they involve considerations of turbulent heat diffusion and turbulent bulk 
viscosity, for which there is currently very little understanding in the context of MRI 
turbulence. One hopes that simulations will shed light on these issues, and currently 
global MRI simulations in the RIAF regime appear to be consistent with thermal 
and viscous stability (see Chapter 2.4). As we discuss further below, simulations of 
the optically thick solutions, which are far more challenging, currently cast doubt on 
alpha prescription stability analyses, even in the supposedly unstable geometrically 
thin, optically thick radiation pressure dominated solution. 

It should also be noted that, in addition to the thermal and viscous instabil- 
ities, other, shorter time scale instabilities exist that are driven by the thermody- 
nam ics of the alpha stress pre scription, particularly the excitation of acoust ic modes 
fe.g. iBlumenthal. Yang fc Linl il984: Kato, Honma & Matsumoto 1988; Chen fc TaamI 
Il993l ). and these might be relevant for explaining high frequency variability in black hole 
sources. Again, however, it is far from clear that the time-dependent thermodynamics 
of the alpha prescription accurately represents the time-dependent thermodynamics of 
MRI turbulence. On the other hand, dyn amical excitation of nonaxisymmetric ac oustic 
waves clearly occurs in MRI turbulence i Heinemann fc Papaloizoull2009all2009br ). 



Each of the three solutions is expected to have a distinct relationship between 
overall radiative luminosity and accretion rate. Geometrically thin, radiatively efficient 
accretion disks are expected to have a luminosity which varies linearly with accre- 
tion rate, as for a given black hole spin, all the released binding energy is equivalent 
to an approximately fixed fraction of the rest mass energy of the accreted material. 
Advective models have reduced radiative efficiency. For the optically thin RIAF solu- 
tions, the luminosity i s approximately proportional to the square of the accretion rate 
(JNaravan fc Yilll995bl ) and therefore the radiative efficiency drops as the accretion rate 
is reduced. Observational tests of this predicted relationship are discussed in Chapters 
3.2 and 5.7. In the high luminosity regime of slim disks, the radiative output is ex- 
pected to approximately saturate to the Eddington limit as photons become trapped 
( Begelman fc Meien 19821). and t his is expected to remain true even if outflows are 



driven (e.g. IPoutanen et al.ll2007h . althoug h much of this luminosity w ill be emitted 
anistropically toward the rotation axis (e.g. lOhsuga fc Mineshigell201ll : Chapter 5.3). 

The three basic flow paradigms have been used over the years to explain the ob- 
served variety of black hole accretion sources. Something like the radiation pressure 
dominated advective slim disks are probably relevant for luminous quasars and QSO's 
("quasi-stellar objects"), narrow line Seyfert I's, ultraluminous X-ray sources, SS433, 
and perhaps some of the intermediate/steep power law states of black hole X-ray bina- 
ries. Geometrically thin, optically thick, radiatively efficient accretion disks extending 
down close to the central black holes are probably most relevant for the high/soft state 
of black hole X-ray binaries and perhaps for some QSO's. Optically thin RIAF's are al- 
most certainly relevant for low luminosity active galactic nuclei (AGN), jet-dominated 
nonthermal AGN such as M87, the Galactic Center source Sgr A*, the inner regions of 
some broad line Seyfert I's, and the inner regions of the low/hard state of black hole 
X-ray binaries. Many of the other chapters in this book address how well these models 
work in explaining the observed properties of these sources. 



2.2 Tying local models into global models: the overall geometry of the flow 



Figure [2] shows that, provided alpha is not too high, the RIAF solutions at any par- 
ticular radius always terminate above a critic al accretion rate where radiative cooling 
becomes 



Ldms always terminate above a critic al accretion rate wnere radiative cooling 
comparable to turbulent dissipation ( Abramowicz et al.lll995l : lNarayan fc Yil 



Il995a ). This critical accretion rate generally decreases with radius at large radii, so 
that provided the accretion rate is not too low, the accretion flow will generally be 
in the radiatively efficient, geometrically thin disk state at large radii. However, at 
smaller radii, provided the accretion rate is not too high (M < fOa Me, where Me 
is the Eddington luminosity divided by c : iNaravan &: Yillf995bl ). the accretion flow 
can exist in one of two distinct thermal equilibria: the optica lly thin RIAF s olution , 
and the radiatively efficient geometrically thin disk solution. iNaravan fc Yil ( 1995bj) 



argue that evaporation from the surface of the thin disk will tend to drive the accre- 
tion flow into the RIAF regime whenever it is possible. One would then be left with 
a geometry which consists of an outer geometrically thin disk extending down toward 
a transition radius inside of which the flow adopts the RIAF solution. As the external 
accretion rate increases, the transition radius moves inward. This then provides an 
explanation for transitions between hard and soft states in black hole X-ray binaries 
( Esin. McClintock fc Naravanlll997l . Figure O. As discussed in Chapter 2.5, it is now 



well-established that the external accretion rate is not the only parameter that controls 
state transitions: hard to soft transitions generally occur at higher accretion rates than 
soft to hard transitions. Nevertheless, this geometry of an outer thin disk and an inner 
RIAF has become a popular model for hard states of black hole X-ray binaries, and of 
certain classes of active galactic nuclei. 

RIAF's are not the only way to produce hard X-rays, however. As illustrated in 
the top most panel of Figure [S] it is possible that a corona containing hot or energetic 
nonthermal electrons exists above and below the geometrically thin disk. This could be 
locally generated by, for example, flares associated with buo yant magnetic field lines in 



a manner analogous to the pro duction of the solar corona ( Galeev. Rosner fc Vaianal 
[l979; Haardt & Maraschi 1991). Alpha disk models in which some fraction of the locally 
generated accretion p o wer is dissipated in an external corona have been developed by 



ISvensson fc Zdziarskil ( I99J). Just as in the sun, the actual geometry of the corona 



could be quite complicated, with multiple coronal patches. 



Starting with the work of lMever. Liu, fc Mever-Hofmeisten ( 2000all2000a) and lRozaiiska fc Czernvl 



(|2000l ). many attempts have also been made to build alpha-based models of outer thin 
disks, inner RIAF's, and coronae that themselves are treated as accreting RIAF's but 
which can sandwich portions of the thin disk. Each of the different flow regions ex- 
change energy and mass through thermal conduction, evaporation and condensation, 
and irradiation. As shown in Figure [4l one can even form inner condensed pieces of 
radi atively efficient thin d isks embedde d inside the RIAF /c orona flow in these models 
f e.g. iMever. Liu, fc Mever-Hofmeister 20071 : iLiu. Done fc Taam 2011 '). While such flow 
geometries may well occur in nature, using the same alpha prescription everywhere, 
especially at high latitudes off the midplane, might be problematic, as discussed briefly 
in section 3.2 below. One hopes that thermodynamically consistent global simulations 
of MRI turbulence may shed light on how transitions between thin disks, coronae, and 
RIAF's actually occur. 
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Fig. 3 Suggested flow geometries involving outer geometrically thin/optically thick disks and 
inner optically thin RIAF's f or var ious observed states in black hole X-ray binaries, from 
lEsin. McClintock fc NaravanI ||1997| ). Here m is the accretion rate scaled by ten times the 
Eddington luminosity divided by c^ . 



2.3 Other variants 



In addition to the three basic local models of optically thin RIAF's, radiatively efficient 
geometrically thin disks, and radiation pressure dominated advective slim disks, nu- 
merous other models have been proposed over the years. Some of these are variations 
on the three basic models, while others involve more significant departures from the 
physics included in these models. 

Because much of the dissipated accretion power is not radiated away, both the 
optically thin and optically thick advection dominated solutions involve fluid which 
is only weakly bound to the black hole, i.e. with internal energy comparable in mag- 
nitude to the orbital binding energy. Outflows are therefore ve ry likely to occur in 
these regimes I Naravan fc Yilll994l : lBlandford fc Begelmanlll999l V and one way of in- 
corporating them (with little physics beyond the invocation of self-similarity) is to 
simply assume that the accretion rate varies as a power law with radius M oc Rf 
(jBlandford fc Begelmanl ll999). Assuming accretion velocities scale with the free-fall 

1/9 

speed oc -R ' , as self-similarity would require, such solutions have a radial density 
proflle p oc BF~ ' , with p = corresponding to the standard ADAF. Such solu- 
tions have been dubbed Adve ction-Dominated Inflow-Outflow Solutions (ADIOS) by 
iBlandford fc BegelmanI ( 19991 ). Outflows have been commonly observed in global sim- 
ulations of these flow regimes, (Chapters 2.4 and 5.3), though exactly how much mass 
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Fig. 4 Possible accretion flow geo metries in evaporation/condensation models, from 
iMever. Liu, fc Mever-Hofmeisteil l l2007f) . 



is lost compared to how much is a c creted remains theoretically uncertain. Very re- 
cent simulations by iNaravan et al.l 1 20121 ) in the low luminosity RIAF regime find 
that outflows are n ot as powerful as previously thought. Time-dependent self-similar 
A DAF solutions by Qgilvid (1999|), that differ from the stationary self-similar solution 



of lNarava 



Yil l|1994f ). appear to be consistent with this result. At high luminosities, 
outflows can also be driven directly by radiation pressure, particularly with the high 
atomic opacities expected for gas around supermassive b lack holes in active galactic 
nuclei (|Murrav et al.ll 19951 : IProea. Stone fc Kallmanll2000l l. 

Another consequence of not dissipating accretion power in advection dominated 
flows is that the entropy of the plasma increases inward, implying that the flow will be, 
at least hydrodynamically, convectively unstable. It has been suggested that large scale 
convection in the flow can transport angular momentum inward, and stationary solu- 
tions in which only small accretion rates occur as material circulates again and again 

1/2 

in convective eddies have been proposed which hav e radial density profiles p oc R ' 
(Convection Dominat ed Accretion Flows or CDAF's: lNaravan. Igumenshchev fc Abramowica 
l200d : iQuataert fc Gru zinov 2000). Th e physical consis tency of such solutions is con- 
tr oyersial (Balbus fc Hawlevll2002l : iNaravan et al]l2002l ). The recent RIAF simulations 
bv lNaravan et al.l (|2012l ) also do not find obvious signs of convection. 

Yet another variant on the optically thin RIAF model is the luminous hot accre- 
tion flow (LHAF) model of lYuanI I 200H ). As we discussed above, unless alpha is large. 



RIAF's have a maximum accretion rate above which the solutions with advective cool- 
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ing do not exist. This creates a problem for using tiiem to explain hard state sources 
that are observed to exist at high luminosities. The LHAF model solves this problem 
by positing that advective cooling is replaced by advective heating. The heating here 
is due to compressional work, and is non-dissipative. LHAF's therefore essentially bal- 
ance compressional work plus turbulent dissipation with radiative cooling. However, 
such models appear to be thermally unstable given their local equilibrium curves on 
the accretion rate vs. surface mass density plane, although the growth rates might be 
long enough to not significantly affect the flow (Yuan 2003 ). 

Models involving strong and/or large-scale magnetic fields in the fiow have also been 
aroposed, f or example flows in which magnetic pressure dominates thermal pres sure 



proposed, t or example nows in wnicn magnetic pressure dominates tnormal pres sure 
(IPariev. Blackman fc Boldvrev 2003: Meier 2005; Machida. Nakamura. & Matsu motd 
I2OO6I : iBeeelman fc Pringld 120071 ) , and flows i n w hich angu lar moment um transport is 



dominated by large-scale MHD outflows fe.g. iFerreira fc Pelletieij|l995h . It is even pos- 
sible for enough magnetic field to be advected inward in the accretion flow that it 
becomes strong enough to disrupt the flow, allowing accretion to only proceed through 
magnetic Ray leigh- Taylor interchange motion s (so-called Magnetically Arrested Disks; 
iNaravan. Igumenshchev fc Abramowica 120031 . see also Chapter 2.4). Such flows are 
beyond the scope of this particular review, which focuses on accretion driven by tur- 
bulence, although they may very well be important in nature. It is now time for us to 
turn to what is known about the local properties of MRI turbulence in accretion flows. 



3 Going beyond alpha: MRI turbulence 



Since the discovery of its rele vance to the physics of accretion flows ( Balbus fc HawlevI 



ll99ll : lHawlev fc Balbuslll99ll ). tremendous theoretical effort has been expended to try 
and understand the properties of MRI turbulence in various regimes of relevance to 
astrophysics. Studying MRI turbulence has allowed us to pose questions that simply 
cannot be asked within the alpha prescription, and has led to considerable new physical 
insight into how the accretion process works, as well as sharpened the true physical 
uncertainties. It is, after all, only by building models based on real physics, rather 
than precriptions that sweep undertainties into a single parameter, that real scientific 
progress can be made. However, it has to be admitted that the ultimate goal of replac- 
ing alpha-based modeling of accretion powered sources with observationally falsifiable 
models based on the actual physics of the turbulence has not yet been achieved. In this 
section I will review some of the fundamental physics issues that have been understood, 
or at least revealed, by local simulations of MRI turbulence. Global MRI simulations 
of black hole accretion flows have also provided considerable insight, and are in fact 
closest to realizing the goal of providing observationally testable models of the accre- 
tion flow onto the Galactic Center black hole source Sgr A*. These global simulations 
are reviewed in Chapter 2.4. 



3.1 Shearing box simulations of MRI turbulence 

It is in the very nature of fully nonlinear, strong turbulence that energy released 
or injected from large spatial scales passes quickly down to microscopic dissipation 
scales through a turbulent cascade. In our case gravitational binding energy is released 
through the MRI which grows by tapping directly into the free energy associated with 
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the differential rotation inherent in the accretion flow. The microscopic dissipation 
scales are associated with the true viscosity and resistivity of the plasma. The actual 
physical dissipation scales relevant to black hole accretion flows are extremely small 
compared to the energy release scales of the MRI (presumably of order the disk thick- 
ness), but numericists have nevertheless hoped that by putting in enough grid zones 
into their simulations, that some convergence can be achieved in describing the prop- 
erties of MRI turbulence. That hope is best achieved in local shearing box simulations 
of the turbulence, where all the computing power is devoted to resolving scales within 
the turbulent cascade, and not on the larger scale dynamics associated with the over- 
all flow geometry (as important as these larger scales are to ultimately understanding 
observed sources). 

The geometry and pro p erties of the shearing box are very nicely described by 
iHawlev. Gammie &: Balbud l|l995l ). Essentially a small, perfectly rectangular Carte- 
sian box is placed in the rotating shear flow, and corotates with the background flow 
at the center of the box. The curvature in the background flow streamlines are entirely 
neglected, but the effects of rotation are nevertheless included through Coriolis forces 
as well as centrifugal forces that are combined with the gravitational force through an 
effective potential. Boundary conditions are such that the flow is assumed to be per- 
fectly periodic in the azimuthal direction, but shearing periodic in the radial direction: 
one imagines many identical shearing boxes sliding past each other according to the 
background differential rotation. If the box is placed in the midplane of the flow, one 
sometimes neglects the vertical gravity and adopts periodic boundary conditions in the 
vertical direction {unstratified shearing boxes), but one can also include vertical gravity 
{stratified shearing boxes) and adopt outflow boundary conditions, or, for computa- 
tional convenience, retain vertical periodic boundary conditions (a stack of accretion 
disk pancakes!). 

The symmetries of the standard shearing box mean that there is no net accretion of 
mass through the box, and therefore in fact no release of gravitational binding energy. 
All the energy associated with the turbulence in fact arises from the net work done by 
the turbulent stresses on the shearing walls of the box. 

Any initial net vertical magnetic field must be conserved in a shearing box sim- 
ulation, and this is also true of net azimuthal magnetic field if the shearing box is 
unstratified (such field can leak out of the vertical boundaries of stratified boxes if out- 
fiow boundary conditions are employed). Shearing box simulations can therefore have 
the net magnetic flux through the box as a flxed external parameter. Unstratified simu- 
lations with no external flux, and no explicit treatment of viscosity and resistivity which 
would resolve dissipation on scales larger than the grid scale, actually produce MRI 
turbulent stresses that monotonically decrease with increa sing numerical resolution 
( Pessah. Chan fc Psaltial2007l : iFromang fc Papaloizoij|2007l )! However, this is a singu- 



lar situation, as including explicit viscosity and resistivity in the MHD equations (albeit 
at far larger values than are relevant for black hole accretion flows) or a net magnetic 
flux does lead to converged levels of stress as numerical r esolution is increased. This 
stress increases with the amount of external magnetic flux ( Hawlev. Gammie fc Balbua 



I995I : IPessah. Chan fc Psaltid [2007|) and also increases with the dimensionless ratio 
of kinematic vis cosity to Ohmic resistivity, known a s the mag netic Prandtl number 



( Lesur fc Longar etti 2007; Fromang ct al. 2007; Si mon fc Haw lev 2009). The resuff 



that stress increases w ith net magnetic flux has also been con firmed in localized regions 
of global simulations ( Sorathia. Revnolds fc Armitagell20ld 'l. 
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Fig. 5 Horizontally-averaged azimuthal component of the magnetic field as a function of 
height z and time in two vert ically stratified shea ring box simulations whose only difference 
is the height of the box, from [Davis, Stone fc Pc ssah (2010). There exists a dynamo in such 
vertically stratified simulations that causes quasipcriodic azimuthal field reversals. 



Adding in vertical gravity in stratified simulations enables convergence of the tur- 
bulent stresses even without net magnetic flux and exp licit viscosity and resistivity 



Duient stresses even witnout net magnetic nux ana exp licit viscosity ana resistivity 
( Shi. Krolik fc Hirosell20ld : IPavis. Stone fc Pessahll201Cl ). A major difference between 



unstratified and stratified shearing box simulations is that gravity allows for mag- 
netic buoyancy, and this is clearly playing a role as alternating signs of azimuthal 
field continually dev elop in a dynamo within the MRI turbulence and rise outward 
( Brandenburg et al.l ,1995'). For weak magnetic fields, this produces a quasiperiodic 
pattern of field reversals, as shown in Figure [S] (When time-reversed, the pattern 
here resembles the butterfiy diagram of latitudunal distributions of sunspots over the 
course of the solar cycle. By analogy, this behavior is occasionally referred to as the 
MRI butterfly diagram.) Moreover, stratifled shearing box simulations with vertical 
outflow boundary conditions and net vertic al magnetic flux can actually locally pro- 
duce magnetocentrifugally driven outflows (ISuzuki fc Inutsukal l2009l : iFromang et al.l 



I2OI2I : iLesur. Ferreira fc Qgilviel I2OI2I : iBai fc Stonel I2OI2I ') . Increasing the net vertical 



magnetic flux in such simulations ca n also destroy the p eriodicity of the MRI butterfly 
dynamo, and ultimately suppress it ( Bai fc Stonell2012l 'l. 



The whole premise of the alpha prescription is that the stresses are inherently lo- 
cal: (vertically-averaged) stress just depends on (vertically-averaged) thermal pressure. 
This appears to be mostly conflrmed by shearing box simulations which have wide 
radial extents, in that spatial correlations between stresses at different locations de- 
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crease rapidly on radial scales larger than the disk scale height. However, there remain 
~ 20 percent correlations in the Maxwell st ress on larger radial scales, indicatin g that 
the turbulence may not be entirely local I Simon. Beckwith fc Arniitagell2012l 'l. The 



butterfly dynamo cycles have also been observed in global simulations and have sig- 



nifica nt radial coherence on scales much larger than a disk scale height ( O'Neill et al.l 

Eonl). 



3.2 Aspects of the vertical structure revealed by MRI simulations 

Provided any external magnetic flux is not too high, stratified shearing boxes generally 
result in a structure that is dominated by thermal p ressure in the midplane and mag- 



result m a structure that is dominated by thermal p ressure m the midplane and mag- 
netic pressure and tension forces in the outer layers ( Stone et al.lll996l : iMiller fc Stond 



l2000l : lHirose. Krolik &: Stonell2006l l . MRI turbulence is generally confined to the weakly 
magnetized region s near the midplane, while P arker instability dynamics dominates 



the outer regions i Blaes. Hirose fc Krolild 120071 ). This basic structure of weakly mag- 



netized midplane regions and more strongl y magnetized high altitude regions is also 
generally observed in global simulation s fe.g. lHawlev fc Balbusll2002l : |Penna et ahlboid : 
ISorathia. Revnolds fc Armitagell201(]f l. although here the strongly magnetized regions 
at high altitude can involve significant radial fiows and circulation, which cannot hap- 
pen in a shearing box. The fact that the magnetic pressure profile is broader than 
the thermal pressure profile, and that Maxwell stresses generally dominate Reynolds 
stresses in the turbulence, implies that an average of the ratio of stress to thermal 
pressure (alpha!) generally increases outwards: alpha should never be treated as a local 
quantity, but instead it is, at best, a representation of the ratio of vertically-averaged 
stress to pressure. Disk atmospheres are generally supported by magnetic fields, not 
thermal pressure, so that if atmosphere models of thermal and reflection spectra rely 
critically on vertical hydrostatic equilibrium between gravity and thermal pressure gra- 
dients, they may not be accurate. In addition, the fact that MRI turbulence is generally 
confined to the weakly magnetized midplane regions suggests that models of accreting 
coronal flows discussed in section 2.2 above are probably not well described by a simple 
alpha prescription. 

The magnetically-dominated outer layers are very suggestive of the locally gen- 
erated magnetized corona discussed above in sec tion 2.2. Indeed, very tall stratified 
shearing box simulations bv lMiller fc Stond (|2000|) found that approximately a quarter 
of the magnetic energy generated in the turbulent midplane regions was carried out 
into the corona. These simulations assumed an isothermal equation of state, however, 
and vertically stratified simulations that capture turbulent dissipation as heat and in- 
corporate diffusive radiation transfer have generally found that the fraction of accretion 
power that is dissipated in the magnetized corona outside the photosphere is very small 
(JHirose. Krolik fc Ston J2006l : iKrolik. Hirose fc Blaesll2007l : iHirose. Krolik fc B laes"2009'') . 



On the other hand, the magnetic buoyancy exhibited in the butterfiy diagram il- 
lustrated in Figure [S] can be very energetically important in transporting significant 
amounts of thermal energy out ward in the fo rm of trapped photons in the radiation 
pressure dominated regime (iBlaes et al.l 1201 il l. It should also be emphasized that the 
existing radiation MHD simulations have not included a net vertical magnetic flux, 
which might in principle enhance the coronal energetics. 
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3.3 Physics issues in the RIAF regime 

Low luminosity RIAF models continue to be plagued by significant uncertainties in 
the microphysics of the plasma, whether they are globally simulated with MRI tur- 
bulence or modeled with an alpha prescription. Because the accreting plasma retains 
a significant fraction of its binding energy as internal energy, temperatures must ap- 
proach virial temperatures: kT ~ GMnip/R = {Rg/R)mpC , where rrip is the proton 
mass. This corresponds to ~ 10 K at ten gravitational radii. Optically thin cooling 
by electrons at such temperatures is very fast, unless the accretion rate and density 
is extremely low, so in order to not radiate away all the heat on an inflow time, the 
electron temperature To must be significantly less than the ion temperature Tj. This 
in turn implies that the electrons should not receive the vast majority of the turbulent 
dissipation of accretion power (otherwise they would radiate it away), and the ions 
must be at least partially thermally decoupled from the electrons on the inflow time. 
Coulomb collisions alone will be insufficient to thermally couple ions and electron s 
provided the accretion rate is not too high ( Rees et al.lll982l : iNaravan fc Yilll995bl ^. 



Plasma instabilities may in principle exist that couple the spec i es more rapidly, but so 
far this has not been demonstrated ( Begelman &: Chiuehl ll988: Par k et al.l 2010l). 



E ven the MRI behaves differently in the collisionless regime (Quat aert. D orland fc HammettI 



I2OO2I : ISharma. Hammett. fc Quataertlliooi : lBalbus|[20o3 : llslam fc Balbudboosi '). espe- 
cially in giving rise to anisotropic pressure te nsors that themselves can give rise to 
significant angular momentum transport (Sha rma et al.l l2006). and this is not cap- 
tured in simulations that assume MHD. How the turbulent dissipation is ultimately 
channeled into heating of the ions and electrons (or energization - the ion and elec- 
tron distribution functions need not be thermal) is another significant uncertainty. 
Local simulations find that direct heating by the anisotropic pressure tensor can ac- 
count for 50 percent of the heatin g by the turbulence , and that the ratio of electron 
to ion heating is ^ 0.3{Tc/Ti) ' ( Sharma et al.l 120071 1. It is also just now becoming 



possible to do fully k inetic simulations of collisionless MRI turbulence, at least locally 
(JRiauclme et al.ll2012 '). so that there is hope for further resolving some of these issues 
in the not so distant future. 

All these effects remain to be included in global simulations of the RIAF regime, 
discussed in Chapter 2.4. Currently these generally assume regular MHD, and adopt 
prescriptions for treating the electron and ion distribution functions, such as assuming 
thermal distributions with a constant ion to electron temperature ratio. 



3.4 Attempts to simulate the radiation pressure thermal instability 

For some years now, it has been possible to do vertically stratified shearing box 
simulations that explore the thermodynamics of MRI turbulence. These simulations 
capture grid scale losses of magnetic and kinetic energy, and incorporate radiation 
transport and cooling through flux limited diffusion. It has also been possible to do 
global simulations using flux limited diffusion under axisymmetry - see Chapter 5.3. 
Recently, even more accurate radiation transport algorithms have been successfully 
developed (e.g. Davis, St one fc Jianeil2012l : |jiang. Stone fc Davisll2012l : ISadowski et al.l 
[2012; Takahashi ct al. 20lli! 

Vertically stratified shearing box simulations with optically thick cooling and which 
incorporate the dynamics of radiation pressure enable the exploration of the thermal 
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instability predicted by alpha disk modeling on the negative slope branch of the ther- 
mal equilibrium curves on the right hand side of Figure [2] The origin of this thermal 
instability is very easy to understand. Assuming that radiation diffusion dominates 
the vertical heat transport (which is generally found in simulations), the local cool- 
ing rate per unit area is of order the radiation energy density at the midplane aT 
times the speed of light over the vertical optical depth. Because the opacity is domi- 
nated by electron scattering in these high temperature regimes, this then implies that 
Q~ oc T /E. On the other hand, the heating rate per unit area is the vertically in- 
tegrated stress times the rate of strain, Q ~ Hw j^^R\dfi / dR\. Because the disk is 
supported vertically by radiation pressure, and the vertical gravity increases linearly 
with height above the midplane, the disk half thickness is simply proportional to the 
surface radiation flux, i.e. H oc Q~ oc T / S. Hence a standard alpha prescription, in 
which Wfl.0 — aP = oiaT /3, will mean that the heating rate Q oc T / E. Because 
the inflow time is much longer than the thermal time for geometrically thin disks, the 
surface mass density E cannot vary significantly on the thermal time scale, and the 
heating rate therefore depends much more sensitively on temperature than the cooling 
rate. Hence a perturbative increase (decrease) in temperature would lead to runaway 
heating (cooling). 

As shown in Figure [ij simulations are consistent with the standard alpha pre- 
scription that the time and space averaged stress scales with total thermal pressure, 
which is mostly radiation pressure in this regime. Nevertheless, at least some sim- 
ulations have been able to establish thermal equilibria between heating and cool- 
ing in the radiation pressure dominated r egime that last for many thermal times 



(ITurnen 2004: Hirose. Krolik fc Blaesll2009h . On the other hand, recent simulations 
bv I Jiang. Stone &: David (2013bj) find that such equilibria, even if established, always 



eventually suffer runaway heating or cooling. It is not clear what is producing the 
differences between the different simulations, which are run using different codes. 

Note that, in contrast to the hydrogen ionization driven thermal/viscous instability 
that is responsible for dwarf nova and outbursts in X-ray binaries, there is no similar 
observational evidence for the putative radiation pressure driven thermal instability 
predicted by alpha disk theory. It could be that the intrinsically stochastic nature of 



the turbulent dissip ation, which is not cap tured in the alpha prescription, acts as a 
stabilizing influence ( Janiuk fc Misrall2012l ). But in addition, it is clear that the alpha 



prescription breaks down on short time scales. As shown in Figure [6] fluctuations in 
thermal energy lag fluctuations in turbulent energy by approximately a thermal time. 
This is easy to understand physically: it is the dissipation of turbulence that heats the 
plasma and produces thermal pressure, and thermal energy will therefore only respond 
to fluctuations in turbulent energy on time scales of order the heating time ~ {afi)~ . 
There appears to be no direct feedback from pressure to stress on the thermal time 
scale, and the alpha prescription is therefore only established on longer time scales, 
although exactly why the prescription is established remains mysterious. More thought 
needs to be applied to understanding the relationship between stress and pressure, 
and to understanding what is going on to e xplain the different results seen in different 
simulation codes. Recently, |Lin_et_alJ (|2012l ') have taken a useful first step in developing 
analytic, time-dependent ordinary differential equations that successfully reproduce a 
number of the observed features in these simulations. 

Note that while the thermal instability has been studied using local simulations of 
MRI turbulence, the inflow or viscous instability has not as the boxes that have been 
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Fig. 6 Cross correlation coefficient for various forms of volume integrated energy as a func- 
tion of time difference with respect to variations in volume integrated magnetic energy, for a 
radiation dominated stratified shearing box simulation. Negative values on the horizontal axis 
mean that the energy lags behind magnetic energy. The dashed line shows the turbulent ki- 
netic energy, which is highly correlated with magnetic energy at zero lag. This is because both 
energies are aspects of the same MRI turbulence! The solid and dotted curves show radiation 
and gas internal energies, respectively. Both of these are very similar as they are both thermal 
energies, and both are again correlated with magnetic energy, but with a significant time lag 
of order 5-15 orbital periods. This is comparable to the thermal time in this simulation. (From 
iHirose, KroHk fc Blae"il200Sl .') 



used are radially too narrow to allow for significant variations in surface mass density. 
Very wide radial boxes, or global simulations, will be necessary to explore this physics. 



3.5 Other issues in the radiation dominated regime 



Radiation dominated plasmas have a number of other interesting properties that are 
likely relevant to black hole accretion flows in the high luminosity regime. One in 
particular is that the sound speed in the fluid ~ [4aT /(9p)] ' is determined by 
radiation pressure, not gas pressure, and when the former exceeds the latter, it is 
possible to be in a regime where turbulent motions are subsonic and yet supersonic 
with respect to the sound speed from gas pressure alone. But photons generally diffuse 
through the plasma, and if they do so rapidly, then even fluid motions that are subsonic 
with respect to the radiation pressure sound speed, but supersonic with respect to the 
gas sound speed, can be highly compressible, because photon diffusion reduces the 
photon pressure response. Large density fluctuations can therefore be produced, and 
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this has been obs erved in unstratified shearinR box simulations of radiation dominated 
MRI turbulence (|Turner et al.ll2003l : [jiang. Stone fc Davisll2013al ). 



Radiation damping of temperature fluctuations in radiation pressure dominated 
MRI turbulence can be a significant source of dissipation which, like the pressure 
anisotropics in coUisionless MRI turbulence discussed above in section 3.3, can be 
resolved in numerical simulations. Such fl uctuations can be compressible in nature 
! just mentioned ( Agol fc Kroliklll998l ). or due to nonlinear isoba ric fluctuations 



associated with regions of high magnetic pressure ( Blaes et al.ll201ll 'l. Some tens of 



percent of the total dissipati on has been observed to occur t hrough radiation damping 
in shearing box simulations l Turneill2004l : iBlaes et al.ll201in . This radiation damping 



can also increase the bulk viscosity and therefore the magnetic Prandtl numb er, and 
can increase the Maxwell stress in the turbulence ( Jiang. Stone fc Davidl2013al ). 



The fact that MRI turbulence can be supersonic with respect to the sound speed 
in the gas alone in the radiation dominated regime implies that, in principle, turbu- 
lent speeds can exceed mean thermal speeds not only of the ions in the plasma, but 
also the electrons. If turbulent motions are limited by the radiation sound speed then 
this may start to happen at radiation to gas pressure ratios in excess of the proton to 
electron mass ratio, and may happen at even lower ratios in the photosphere regions 
which tend to be dominated by magnetic pressure, not thermal pressure. Differences in 
bulk turbulent velocities on the scale of a photon mean free path that exceed in mag- 
nitude characteristic electron thermal speeds will mean that bulk Comptonization by 
the turbulence itself will dominate thermal Comptonization, and this may provide an 
alternative means of producing a Comptonized high energy spectrum in radia t ion dom- 
inate d luminous states of black hole sources ( Socrates. Davis fc Blaeal2004l : ISocrated 

boiol i. 



Advection of heat is a key ingredient to the radiation pressure dominated slim disk 
solutions discussed in section 2.1 above. However, instabilities might produce inhomo- 
geneities in the flow that allow photons to escape more readily through underdense 
channels, rather than be advected inward. The most well- explored of s uch instabilitie s 
are magnetically-mediated "photon bubble" instabilities (: Aronal 19921 : lGammielll998l l. 
On short length scales where photons diffuse rapidly, such inst abilities amount to radia- 
tively amplified magnetosonic m odes (IBlaes fc Socrateal2003l') that d evelop into highly 



inhomogeneous trains of shocks ( Begelmanll200ll : lTurner et al.ll2005l . see Figure [7|. In 



principle such inhomogeneities could allow locally supe r-Eddington flux es to escape 
from the disk atmosphere without driving an outflow ( BegelmanI [2003) • Simulating 



photon bubbles in the presence of MRI turbulence has proved computationally chal- 
lenging in the radiation pressure dominated regime due partly to the small length scales 
(of order the gas pressure scale height) that must be resolved, and partly by the fact 
that Parker instabilities in the magnetically dominated surface layers also produce sig- 
nificant inhomogeneity. Mo dels of slim disks with p orous outer layers and winds have 
recently been developed bv iDotan fc Shavivl (|201lh . It has also been suggested that 
accretion flows in the radiation dominated regime might be highly inhomogeneous 
structures that are not well-described by an y of the standard accretion flow models 
discussed in section 2.1 ( Dexter fc Agolll201ll ). 
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Fig. 7 2D simulations of shock trains produced by the photon bubble instability in radiation 
pressure supported media with initially uniform magnetic fields of increasing strength from 
left to right. Arrows show fluid velocity and colors show density on a logarithmic scale, with 
warm colors being high and cold colors being low. The weight of the dense shocked fluid causes 
the weaker magnetic flelds in the left hand flgure to buckle. fFrom iTurner et al.ii2005i .) 



4 Summary 



The primary reason why the alpha prescription continues to be a mainstay of black 
hole accretion theory is that it enables models to be built that couple the dynamics 
of the flow (outward angular momentum transport through the plasma by a stress de- 
scribed phenomenologically by alpha) to the thermodynamics (dissipation of accretion 
power described phenomenologically by the same alpha stress prescription times the 
shear rate). Once this is all combined with radiative cooling processes, one has models 
that can be used to generate spectra, time variability, etc. that can be compared to 
observations. While decades of theoretical work within this research paradigm have 
yielded valuable insights (e.g. scalings of luminosity and temperature with accretion 
rate, the importance of advection of heat), one can only carry this program so far 
without addressing the fundamental physics that the alpha prescription hides. 

(Fewer) decades of research have now been spent on understanding the properties 
of MRI turbulence, but until recently, this has focused mostly on the dynamics of the 
turbulence, not the thermodynamics, and it is the latter which is ultimately required 
to connect to observations, which are, after all, detecting the photons emitted by the 
source. We still have a lot of unanswered fundamental questions. For example, what 
is the true nature of the thin disk/RIAF transition radius, upon which hangs so much 
phenomenology of black hole X-ray binaries (hard/soft state transitions, band-limited 
noise, low frequency quasi-periodic oscillations, ...)? The thermal/viscous instabilities 
have guided theoretical effort to exclude uphysical equilibria, but are these instabilities 
always real, and how do they manifest? (Only those driven by ionization/recombination 
are observationally known to exist, e.g., in dwarf novae.) What determines how accre- 
tion power is partitioned into various forms? These are all questions of thermodynamics, 
not just dynamics. 

As discussed in Chapter 2.4, global MRI simulations of low luminosity RIAF's have 
been most successful in connecting to observations, largely because, at least in some 
regimes, the radiative cooling does not dramatically affect the dynamics of the flow, and 
therefore can be calculated after the fact by post-processing the simulation. Moreover, 
because the flow is optically thin at most photon frequencies, the simulation hardware 
has actually developed to the point where radiative coolin g can be fully incorporated 
in the actual dynamical simulation itself i Dibi et al.l 120121 ). However, as discussed in 
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section 3.3 above, even here there remain issues in the microphysics that still need to 
be understood theoretically. 

Optically thick radiation transport is computationally more expensive, and also 
more important for the dynamics in the high luminosity black hole accretion regimes 
where radiation pressure plays a critical role. As discussed in section 3.4 above, lo- 
cal shearing box simulations have shed some light on how MRI turbulence works in 
this regime, and in what ways the alpha prescription does, and does not, describe the 
physics. Unfortunately, this still leaves major unanswered questions as to the global 
structure of the flow. Two dimensional (axisymmetric) global simulati ons have also 
been done whi ch have confirmed the existence of discrete flow states I Qhsuga et al.l 



1200 9: Qhsuga fc Mineshigell201ll . see also Chapters 2.4 and 5.3), but these cannot be 
run over long time scales as MRI turbulence cannot be sustained in axisymmetry. 
With the ongoing increase of computer power, combined with the development of new 
radiation transport algorithms, it should be possible to do global 3D simulations of 
accretion flows in optically thick regimes too. This will better enable us to understand 
the origins of state transitions in black hole X-ray binaries, transitions between radia- 
tively efficient thin disks and RIAF's, and other fundamental problems in black hole 
accretion flows. Achieving this goal should complete the shift to a more physics-based 
research paradigm. 
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